


function y=b1(x);

r=0.03;

y=0.141./( r*cutoff1(x) ); 
%y=0.141480./( r*cutoff1(x) );

%It seems I used the rounded version for the constant cutoff case and the
%unrounded version for the variable cutoff case

